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Abstract 

Following Sanchez's approach we investigate the effect of scalar mass in the 
absorption and emission problems of Ad Schwarzschild black hole. The ab- 
sorption cross sections for arbitrary angular momentum of the scalar field are 
computed numerically in the full range of energy by making use of the ana- 
lytic near-horizon and asymptotic solutions and their analytic continuations. 
The scalar mass makes an interesting effect in the low-energy absorption cross 
section for S-wave. Unlike the massless case, the cross section decreases with 
increasing energy in the extremely low-energy regime. As a result the univer- 
sality, i.e. low-energy cross section for S-wave is equal to the horizon area, 
is broken in the presence of mass. If the scalar mass is larger than a critical 
mass, the absorption cross section becomes monotonically decreasing func- 
tion in the entire range of energy. The Hawking emission is also calculated 
numerically. It turns out that the Planck factor generally suppresses the con- 
tribution of higher partial waves except S-wave. The scalar mass in general 
tends to reduce the emission rate. 
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I. INTRODUCTION 



Black hole is an extremely important physical system in which we can test our formulation 
of quantum gravity and string theories. For example, recent understanding of Bekenstein- 
Hawking entropy [1,2] based on the microstates [3,4] gives an impetus to recent drastic 
development of the non-perturbative string theories. The typical phenomenon appearing in 
black hole physics as a quantum gravity effect is an Hawking radiation [5,6], which states that 
black hole emits a thermal radiation at Hawking temperature. The Hawking temperature 
T BH for the non-rotating and chargeless black hole is T BH = l/A-nr s where r s is an horizon 
radius. For example, the emission rate of energy in each mode of frequency k and angular 
momentum i can be written as 



where the upper (lower) sign in the denominator is for boson (fermion). The quantity T^{k) 
in numerator is an absorption coefficient, which is also named 'greybody' factor. From a 
viewpoint of quantum mechanical scattering theory, it is related to the scattering amplitude 
Se(k) in the form 



In usual quantum mechanical systems Si(k) is an unitary quantity, i.e. Si = e 2lSe , where 
the real parameter 5i is a phase shift. Thus, T^k) vanishes in these systems. However, the 
amplitude Se(k) is not in general phase factor in the black hole system [7], which makes the 
phase shift to be complex, i.e. 5e = + ify. The non-unitarity of the scattering amplitude 
might be closely related to the information loss problem of the black hole [8-10] although 
Hawking stated the information loss puzzle in terms of the quantum pure state and quantum 
mixed state [8]. 

Eq.(l.l) indicates that the calculation of the absorption coefficient T t (k) is essential 
for the computation of the Hawking radiation. As will be shown shortly (see Eq.(2.5) and 
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Eq.(4.3)), Te(k) can be straightforwardly computed from the partial absorption cross section 
<Te(k). Thus, it is important to compute at(k) strictly as much as possible. 

The rough estimation of the absorption and emission rates for the various black holes 1 
were computed long ago [12-15]. Especially Unruh derived in his seminal paper [16] the 
analytic expression of the low-energy partial absorption cross section for the massive scalar 
and Dirac particles by matching the near-horizon radial solution and the asymptotic solu- 
tion via the solution in the intermediate regime. It is instructive to introduce an explicit 
expression for the partial absorption cross section of Ref. [16] for the massive scalar: 

y^i) unruh "TT 

^{2£ + l)T e (1.3) 



ir{£\) 4 2 2e+2 {l+v 2 )k 2e+3 v 2e ^ 



(2£!) 4 (2£ + 1) 2 [1 - exp {-7rfc(l + v 2 )/v}] 
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' k(l + v 2 ) \ 2 
2v ) 



where v = yl — m 2 /k 2 . When Unruh derived Eq.(1.3), he put the horizon radius of the 
Schwarzschild black hole to be unity. 

The most interesting result in Eq.(1.3) is the S-wave cross section which is 

(rr ) - (27T) 2 (l+* 2 )fc 

Taking m — * limit and assuming k << 1, we have (<7o) unru h ~ 47r, which equals to the 
horizon area. This is an universal property of the low-energy absorption cross section for 
massless scalar in black hole physics. 

Inspired by string theories, the universal property is re-examined few years ago in the 
arbitrary-dimensional black hole [17] and higher-dimensional objects [18] such as p-brane. 
More recently, it is also examined whether or not the effect of mass maintains the universal 
property of the low-energy absorption cross section in the higher-dimensional theories [19,20]. 
Although the mass effect in general does not break the universal property on condition that 



The Hawking radiation for near-extremal 5d black hole and its implications in the information 
puzzle are discussed in Ref. [11] 
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mass is not too large, it is found that the mass-dependence of the cross section is very- 
sensitive to the spacetime property especially, in the near-horizon regime. 

While the low-energy behavior of the absorption and emission in Ad black hole and higher- 
dimensional objects were extensively examined, the behavior in different range of energy is 
not sufficiently examined. Fortunately, for _D3-brane system in Born-Infeld theories [21,22], 
however, the wave equation is represented in terms of the modified Mathieu equation. Using 
the various properties of the equation it is possible to examine the absorption cross section 
in high-energy domain as well as low-energy domain [23-27]. However, fairly speaking, it 
seems to be far away to understand the absorption and emission spectra in the full range of 
energy in the general higher-dimensional systems. 

For Ad black hole system the full absorption spectrum of Schwarzschild black hole for 
the massless scalar was computed long ago [28,29] by making use of the analytical series ex- 
pressions of near-horizon and asymptotic solutions and their analytic continuations. In fact, 
these expressions and their mathematical properties were examined in detail by Persides 
[30-32]. In Ref. [29] Sanchez matched these two analytical solutions via analytic continu- 
ation, which enables to compute the absorption and emission spectra in the full range of 
energy. 

In this paper we would like to study the effect of mass in the absorption and emission 
spectra of black hole in the full range of energy. This seems to be important due to the fact 
that most matters have their own masses. However, the effect of mass can be negligible in 
the realistic situation because matter's mass is in general extremely small compared with 
black hole's mass. Our particular interest is to check whether the mass effect maintains 
the universal property of low-energy absorption cross section or not. Also, it is of interest 
to check how the mass effect changes the global properties of the absorption and emission 
spectra. 

The numerical calculation shows that the absorption cross section for massive scalar case 
is in general larger than that for massless case. In particular, for the S-wave the introduction 
of mass makes an unexpected effect in the low-energy regime of the partial cross section. In 
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the extremely low-energy regime the partial cross section decreases with increasing energy 
unlike the case of massless scalar field. This indicates that the universality is generally broken 
due to the scalar mass. As Sanchez found in the massless case [28,29], the total absorption 
cross section has an wiggly characteristic which implies that each partial absorption cross 
section has a peak in different energy scale. However, the total emission rate does not have 
the oscillatory pattern, which indicates that the Planck factor in general suppresses the 
contribution of the higher partial waves. It turns out that the scalar mass generally reduces 
the emission rate. 

The paper is organized as follows. In section II we will describe the results for the 
massless case briefly without any calculation. In the same section we will comment the 
possible physical mechanism which explains the results roughly. In section III the wave 
equation for the massive scalar is solved to derive the analytic solutions in the near-horizon 
and asymptotic regimes. These two solutions are matched with each other by making use of 
the analytic continuation in section IV. The numerical calculation is performed in the same 
section to examine the effect of the scalar mass in the emission and absorption spectra of 
black hole. In section V a brief conclusion is given. In appendix the Wronskians for the 
analytic solutions are explicitly calculated. 



II. BRIEF DESCRIPTION OF RESULTS FOR MASSLESS CASE 



In this section we will comment briefly the results of the absorption and emission prob- 
lems for the massless scalar case. The minimally coupled scalar field equation □$ = in 
the background of Schwarzschild spacetime, whose metric is 



= diag 



-i 



i 2 • 2 a 
,—r ,—r sin U 



(2.1) 



where r s is an horizon radius, reduces to the following radial wave equation 



, cPR, 



dR f 



x(x - x s ) —— + (x - x s )(2x - x s )— h 



dx 2 



dx 
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+ l)(x - X s ) 
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where $ = Re(r)Ye >m (9, (f))e~' lkt , x = kr, and x s = kr s . This radial equation was analytically 
solved in Ref. [30] at near-horizon and asymptotic regions separately as a series form and 
the mathematical properties of them were analyzed in Ref. [31,32]. These two solutions 
are matched in Ref. [28] by the analytic continuation which makes a common domain of 
convergence for the solutions. As a result, one can compute the absorption and emission 
spectra numerically. We do not want to repeat the calculation here. What we want to do in 
this section is just a brief description of the results and possible mechanisms for explaining 
why these results come out. 

Fig. 1 shows the ^-dependence of the partial absorption cross section ai for £ — 0, 1, 2 
in the full range of k when r s — 1. The global shape of <j£ exhibits an increasing and 
decreasing patterns at small k and large k respectively, and each has a unique global 
maximum at k = k\. This can be understood as following. If we use a "tortoise" coordinate 

= x + x s \ia(x/x s — 1), the radial equation (2.2) is changed into the usual Schrodinger 
form as following: 

~^ + K ( />*)^ = ^ (2.3) 
where tp£ = xRg, = kr\ and 




Thus, V}ff and k 2 play roles of quantum mechanical potential and energy eigenvalue respec- 
tively. The potential V^j is plotted in Fig. 2 for I — 0, 1, and 2 with respect to r*. The 
overall shape of is a potential barrier which separates the near-horizon and asymptotic 
regions, and the barrier height is denoted by ujJ. The numerical result of k\ and ue are given 
at Table I. 



Table I 





£ = 


£= 1 


£ = 2 


£ = 3 


£ = A 


k\ 


0.239 


0.687 


1.099 


1.503 


1.903 




0.325 


0.630 


0.994 


1.370 


1.750 



6 



The Table shows kg > oug for nonzero i. This can be understood from the fact that og 
is proportional to the transmission coefficient (or absorption coefficient) and the maximum 
transmission coefficient generally occurs when energy is higher than potential barrier height. 
This fact also roughly explains the increasing behavior of og in low-energy region. The 
decreasing behavior of og in high-energy region can be explained from the relation of ag{k) 
with Yg{k). For massless case this relation is given by 

ag(k)^^(2£ + l)rg(k). (2.5) 

Since Tg(k) goes to unity at most in high-energy region, we have <Je(k) ~ l/k 2 , which explains 
the decreasing behavior of ag in the same region. The fact lj > k$ may occur because of 
the competition between k 2 in denominator and T (k) in numerator of Eq.(2.5). 

Fig. 1 also indicates a (k = 0) = 4n, which is an universal property of low-energy 
absorption cross section for S-wave. The fact of a (k = 0) ^ and agik — 0) = for 
nonzero I implies that r (A:) ~ k 2 and T e (k) ~ k a with a > 2 at low-energy region. This is 
in good agreement with Starobinsky's formula [12,28] 

r '~(2WWTTW {2kr - )2t+2 (2 ' 6) 

at \ow-k regime. The results of Unruh, (ag) unr uh in Eq.(1.3) is plotted by red lines in Fig. 
1. Thus the cross sections {ag) unru h derived by Unruh in Ref. [16] are correct only in the 
extremely low-energy regime. 

The black line in Fig. 3 is a plot of /c-dependence of total absorption cross section, i.e. 
&tot = J2e a e- The wiggly pattern of a to t implies the wave-black hole interaction in the non- 
hermitian effective potential [7]. This also indicates that each partial cross section has a 
peak in different value of k. Similar oscillatory behavior also appears in the total absorption 
cross section of D3-brane [26]. The limiting value at large-A; region is 277r/4 ~ 21.2. Based 
on this limiting value Sanchez made an formula for a tot : 

27n / -sin(v / 27r7r/c) , n „. 

<rtot(k) = — - V2^^-. '-. (2.7) 
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The black line in Fig. 4 is a plot of k- dependence of emission rate, which is computed 

by 

dk [ ' 7T 2 [e 4 ^ - 1] ' 

The oscillatory pattern of a tot disappears in the emission rate because the S-wave contri- 
bution predominates in Hawking radiation. The rapidly decrease of the Planck factor for 
k > 1 suppresses the contribution of the higher partial waves. 

III. ANALYTIC SOLUTIONS FOR THE MASSIVE SCALAR CASE 

In this section we will derive the analytic solution for massive scalar field coupled to 
the Schwarzschild background in the near-horizon and asymptotic regimes. The minimally 
coupled field equation (□ + m 2 )$ = in the Schwarzschild spacetime (2.1) provides a 
following radial equation 

+ nd 2 Ri , , , dRg , . 

x(x - x s ) —— + {x- x s ){2x - x s )—— (3.1) 



+ 



dx 2 dx 

,2 



X 3 - t{t + 1){X- X 8 ) + j^jXsX 2 



Re = 



where v = yjl — m 2 /k 2 , x = kvr and x s = kvr s . Of course, Eq.(3.1) simply reduces to 
Eq.(2.2) when m — 0. Unlike Eq.(2.2), however, Eq.(3.1) involves an explicit m-dependence, 
which makes a situation more complicate. 

At this stage it is worthwhile noting the meaning of r s = 1, which we use frequently 
through a paper. Since the Compton wavelength is inversely proportional to the mass, 
m ~ 1 in this unit corresponds to such a mass of the field when its Compton wavelength is 
order of the gravitational radius. 

Before solving Eq.(3.1) it is instructive to consider the effective potential. Letting x* = 
x + x s \a.{x/x s — 1) and ip£ = xRg make Eq.(3.1) to be 
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d 2 ^e , T/ , ,22, 

+ Veff,lVt = k v ip t 



(3.2) 



V effA r *) 



1 - 



m 2 r. 



where r* = x*/ fci>. Thus V e ffj and /c 2 w 2 play same roles with a quantum mechanical potential 
and energy eigenvalue. 

The ^-dependence of V e ff,t is plotted in Fig. 5 for £ = 0, 1 and 2. The mass term in 
Veff,t generally lowers the barrier height. Thus, there exists an critical mass m c for each 
£ such that the barrier height is lower than the asymptotic value of V e ffjL- The numerical 
calculation shows m c = 0.385 for £ = 0, m c = 0.795 for i — 1, and m c = 1.277 for £ = 2 
when r s = 1. Thus if m > m c for each £, the absorption coefficient T e (k) becomes large, 
which should enhance the absorption cross section in the low-energy domain k rs to. Since, 
furthermore, the mass term in V e f$,t makes the near- horizon regime to be more stable than 
the asymptotic region, Te(k) may become non-zero at k = m ^ 0, which also enhances the 
absorption cross section in the low-energy domain. 

Following Ref. [30] it is straight forward to derive the analytic solutions Qe(x,x s ) in the 
near-horizon region and J-i{±){x,x s ) in the asymptotic region as a series form: 



g t (x,x.) = e -*¥M*-*.l £ den ( x _ Xs y 

n=0 

T m {x,X s ) = (ii^l^x+Vlnl*-*. I) ^ Tn(±)X -(n + l). 



(3.3) 



ra=0 



The coefficients c^ )n and r n (±) satisfy the following recurrence relations: 



n — % — +x.[l + 



to 



k 2 v 2 



x s de,n + 



■ X s 

n — % 1 

v 



.X, 



n 



(3.4) 



+x 3 + 



2m 2 
k 2 v 2 



de, n -\ + 3 



to 

k 2 v 2 



X s de,n-2 + G^,n- 







±2m + x s ^- - 1^ 


Tn(±) - 





(£ + n)(£-n + l)^^ --1 



T"(n-1)(±) 

-(n - l) 2 a; s r (n _2)(±) = 

where To(±) = 1. We will show shortly that the coefficient c^o is an important quantity 
for the study of absorption problem from the viewpoint of quantum mechanical scattering 



theory. In Eq.(3.3) we choosed only ingoing wave in near-horizon solution and ingoing and 
outcoming waves in asymptotic solution. It is interesting to note J~£i±\ = J~ (.{+)■> which comes 
from the real nature of Eq.(3.1). In the appendix we will show the Wronskians become 

2i\ gi \ 2 x 2 s 



W [G}{x, x s ), G e (x, x s )] x = G\d x Qi - Gid x G\ 



vx(x — x s ) 



(3.5) 



W 



^e(+)(x,X s ),J 7 e r_ ) (x,X s 



2i 



-)O x J-t{+) 



x(x — x s ) 
where g t = d efi . 

Now, we would like to show how the coefficient gt is related to the partial scattering 
amplitude. From the viewpoint of quantum mechanical scattering theory the real scattering 
solution Rt(x,x s ) of Eq.(3.1) should be only ingoing wave at the near-horizon region and 
mixture of ingoing and outcoming waves at the asymptotic region as following: 



7 t X >X S 

tie ~ 



(X - X s 



x—>oo -P-l-l 

Ro ~ i + 



21 + 1 
2x 



^[l + 0(x-x s )\ 



(3.6) 



+ - 



X/ 



where Si is a partial scattering amplitude. If we define a phase shift 5g as St = e 2lSe , the 
second equation of Eq.(3.6) can be written as 



Rt 



21+1 



e iSe sin 



x 



x + — In \x — x s \ £ + 5e 

v 2 



+ 



x/ 



(3.7) 



Next let us consider the Wronskian W = R\d x Ri — Rtd x R\. Although we do not know 
the exact expression of Re in the full range of x, it is enough to use the near-horizon 
solution given in Eq.(3.6) for the computation of W, which yields W = —2i\ge\ 2 x 2 /vx(x — 
x s ). Thus same value should be obtained if we use the asymptotic solution (3.7) for the 
computation of W. If, however, the phase shift 5e is real, W becomes zero at the leading 
order which makes a contradiction. Thus we should assume that 5e is complex, i.e. 5e = 
r\i + i@e- Then the Wronskian W computed by the asymptotic solution becomes W = 
—i(2i + l) 2 sinh 2/9£e" 2/3f /x(x — x s ). Therefore, equating those Wronskians gives a relation 



\9t\ 



1 \ 2 

t + h) ( 



-4f3 e 



(3.8) 
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As a result we can obtain gi from the partial scattering amplitude or vice versa. 

Now, we would like to discuss how ge can be computed by matching the asymptotic and 
near-horizon solutions with each other. For this it is convenient [28] to introduce a new wave 
solution <fi(x,x s ), which differs from Re(x,x s ) in its normalization. It is normalized in such 
a way that 

<p t (x, x s ) x ~ Xs (x - x s )-^ [1 + 0(x - x s )\ . (3.9) 

Since J 7 e(±)(x , x s ) in Eq.(3.3) are two linearly independent solutions of the radial equation 
(3.1), one may write tp£ as a linear combination of them: 

ipe(x,x s ) = f^~\x s )J 7 ew (x,x s ) + f\ + \x s )T^-){x,x s ) (3.10) 

where the coefficients /^(a; s ) are called the jost functions. Eq.(3.5) enables us to compute 
fj^ if we know (ft as following: 

(±) ^ xjx - x 9 ) Jxr r<r -I ^ kvrjr -r s ) { i 
it = ± 2i L * (±) ' Vt ix = 2i ^ (±) ' ^Jr ( ^ 

where W[/i,/2] r = f\d r f2 — jid r j\- Inserting the explicit form of ^(±) given in Eq.(3.3) 
and comparing it with the asymptotic expression of Ri in Eq.(3.6), one can easily derive the 
following two relations 



Se(x s ) = J -j^ (3.12) 



and 



Inserting (3.13) into (3.8), therefore, yields 

„2 



9i(x s ) = „(_) ~\ - (3.13) 



1 - \S e Xx s )\" ee 1 - e-^ = * (3.14) 

V\fl ( X s)\ 

The partial absorption cross section is defined in terms of ge(x s ) as [7] 

(3.15) 

21 + 1 y J 
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Combining (3.8), (3.13) and (3.15) we can derive 



1) (l - e~^) = . (3.16) 



TTV 

° l = ^ ' |/;-'i 2 



Thus one can compute all physical quantities such as fy, St and at from the jost functions 

f (±) 
Je ■ 



IV. NUMERICAL CALCULATION 

In this section we would like to discuss the effect of scalar mass in the absorption and emis- 
sion problems. For this we should compute the jost functions /^(x s ) defined in Eq.(3.10) 
numerically. Once f^\x s ) are computed, the partial absorption cross section at is straight- 
forwardly computed by Eq.(3.16). 

The computational procedure is as following. Firstly we should note that <pt(x,x s ) de- 
fined in Eq.(3.10) can be expressed as Qe(x,x s ) in Eq.(3.3) if we assume de t o = 1, i-e. 
ipt(x,x s ) = Gt(x, x s )\d e =i- This is a convergent expansion around the near-horizon region. 
Of course, the other expression of ipt(x,x s ) in Eq.(3.10) is a convergent expansion around 
the asymptotic region. Since, however, the domains of convergence for these two different 
expressions are different from each other, we can not use them directly for the computa- 
tion of /^(xs). In other words, we need two expressions which have common domain of 
convergence. This is achieved by analytic continuation. 

In other to perform the analytic continuation we need a solution of the radial equation 
(3.1), which is a power-series expansion in the neighborhood of an arbitrary point x = b. 
Indeed such a solution can be obtained as following: 

oo 

R(x) = e -*^H*-*«l £ c n {x - b) n (4.1) 

n=0 

where C n satisfies the recurrence relation 



by\n + l)(n + 2)C n+2 + y(n + 1) 



+ 



y(n + 1) + b ( 2n + 1 - 2i— 

\ v 

(b + 2 y)n 2 + jy - 2i^(y + 6)} n + b (b 2 - £\ - y + 



C n+l (4.2) 



x s 1 m 2 b 2 x<, 
1 + i — 



v J k 2 v 2 
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+ 



[ n -l)(n- 2i^j + 3b 2 - £(£ + 1) + (i^- - l) 



2m 2 bx s 

-1 + 



k 2 v 2 



Cn-l 



and y = b — x s . Thus one can compute in principle all the C n 's (n > 2) in terms of Co 
and Ci, which are also expressed in terms of i?(6) and d x R(b). This analytic continuation 
should be repeated from ipi(x,x s ) = Ge(x,x s )\d eo =i to raise the domain of convergence by 
using (4.1) and (4.2). Of course, one can apply a similar procedure to (pi(x,x s ) in Eq.(3.10) 
to lower the domain of convergence. We made a computer program such that the number of 
analytic continuation N is dependent on k. For example in the low-energy region (k ~ to), 
N becomes large because the domains of convergence for each solution are too far. Once we 
make two solutions which have common domains of convergence, the jost functions can 
be computed by Eq.(3.11). 

Fig. 6 show the /c-dependence of fc> r £ = 0,1,2 and various to when r s = l 2 . As 
in the case of massless scalar fa is a monotonically increasing function with respect to k 
for the massive case. However, the scalar mass generally decreases the increasing rate of 
fa. Since fa becomes smaller when to becomes larger for fixed k and £, the factor v/x 2 s in 
Eq.(3.16) increases and 1 — e~ 4A? decreases. Thus, the partial cross section an is determined 
by competition of these two factors. 

Fig. 7 shows the numerical results of /c-dependence of a e for £ = 0, 1, and 2 when r s — 1. 
Fig. 7a shows that the scalar mass tends to increase the absorption cross section for S-wave. 
This means the factor v/x 2 s in Eq.(3.16) predominates the factor 1 — e _4ft in the entire range 
of k. Fig. 7b and Fig. 7c show that the factor 1 — e~ 4ft is superior in the low-energy regime 
for nonzero £ if the scalar mass is not too large. 



2 The unit r s = 1 , which we used in the plots of most figures, implies that the black hole's mass 
is 1/2 in the Planck unit. Thus, the comparability of the field's mass to the Planckian mass in the 
figures represents pure theoretical interest. 
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The most interesting behavior appears at Fig. 7a which is a plot of c^ =0 . For to = 
the low-energy absorption cross section exhibits the manifest universality, i.e. cr ~ 4n. 
However for the massive case this universality seems to be broken. The low-energy region is 
distinguished by two regions on the condition of to < m c . In the extremely low-energy region 
the absorption cross section a is rapidly decreasing with increase of k. In the remaining 
region <j increases with increase of k as usual massless case. As a result, ae has a local 
minimum in the low-energy regime. When, however, the scalar mass to approaches to the 
critical mass m c ~ 0.385, this kind of distinction in the low-energy regime disappears and 
the absorption cross section becomes monotonically decreasing function in the full range of 
k. This seems to be obvious because the potential barrier, which separates the near-horizon 
region from the asymptotic region, is any more meaningless for to > m c . The effect of scalar 
mass in a\ and a 2 is shown in Fig. 7b and Fig. 7c. The cross section with t ^ seems to 
exhibit a similar behavior when the scalar mass approaches to the critical mass. However, 
no distinction in the low-energy regime seems to occur even if to <C m c . Even in this case 
the absorption cross section is increasing function in the low-energy regime with respect to 
k. 

Fig. 8 shows the mass-dependence of S-wave absorption cross section in the low-energy 
limit (k ~ to). In the real calculation we choosed k = to+0.004 and r s = 1. This figure shows 
how the universality of the low-energy absorption cross section is broken due to the effect 
of scalar mass. Fig. 8 indicates that the universality is roughly preserved when to < 0.02. 
However the large mass manifestly breaks the universality, which is also consistent with Fig. 
7a. The red line in Fig. 8 is a plot of TO-dependence of (cr ) unru h, which coincides with <r in 
the neighborhood of the massless region. 

Fig. 3 is a plot of total absorption cross section a tot for various to. In the low-energy 
region the total absorption cross section increases with increase of to, which is also evident 
from Fig. 7. The divergence of the total cross section in the low-energy regime can be 
understood easily from Eq.(1.4). Of course, the massless limit (v — 1 and k ps 0) of (co) unruh 
is 4?t as expected. In the massive case, however, the low-energy limit (v ~ and k to) of 
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it is 4ir 2 m/v 2 , which goes to infinity. Thus, total cross section should go to infinity at this 
regime. In high-energy region a tot becomes m-independent and approaches to the massless 
limit 27tt / 4. This fact reflects an approximate conformal symmetry in the high-energy limit 
[33]. 

Finally we would like to discuss the effect of scalar mass in the emission problem. For 
the massive scalar case the relation between the partial absorption cross section ae and the 
absorption coefficient T e is given by [16] 

Thus the difference from the massless case is a factor v 2 in the numerator. Combining 
Eq.(l.l) and (4.3) makes the emission rate to be 



dH e k 3 v 2 r 2 <7£ 



(4.4) 



dk vr 2 (e 4 ^-l)' 

The total emission rate times 10 4 is plotted in Fig. 4 for various m. As in the massless 
case the oscillatory pattern disappears in the emission rate because of dominance of the 
S-wave contribution. It is interesting to note that in contrast to the absorption case the 
emission rate decreases with increase of m because the factor v 2 in Eq.(4.4) reduces the 
partial emission rate. The Planck factor seems to suppress the contribution of higher partial 
wave even in the masive case. 



V. CONCLUSION 

The effect of scalar field mass is examined in the absorption and emission problems of 
the Ad Schwarzschild black hole in the entire range of energy. The most interesting feature 
may occur in the low-energy absorption cross section for S-wave (see Fig. 7a). The effect 
of mass seems to divide the low-energy region into two parts. In the extremely low-energy 
region (region I) the S-wave absorption cross section er rapidly decreases with increase of 
energy. The decreasing rate becomes larger and larger as the scalar field becomes heavier 
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and heavier. In the moderately low-energy region (region II) <t increases with increase of 
energy as usual massless case. This type of behavior does not appear in massless case, 
whose o"o has only region II in the low-energy regime. When the scalar mass is less than the 
critical mass m c ~ 0.385, where the barrier height of the effective potential (3.2) becomes 
zero, both regions co-exist, but region I extends and region II shrinks with increasing m. 
When m approaches m c , region II finally disappears and the S-wave cross section becomes 
the monotonically decreasing function in the full-range of energy. 

Another interesting feature in the absorption problem is a fact that the heavy scalar 
particle generally breaks the universal property, i.e. lim^o°o ~ S, where S is an area of 
the horizon. The numerical calculation indicates that this universality is roughly preserved 
when m < 0.02 (see Fig. 8). 

There appears a wiggly pattern in the total absorption cross section as in the massless 
case [28,29], which indicates a wave-black hole interaction (see Fig. 3). This oscillatory 
pattern also indicates that each partial absorption cross section has a peak in the different 
point of energy. Athough the total absorption cross section for massive case is larger than 
that for massless case in the low-energy region, both approach to 27n/A in the high-energy 
limit, which reflects an approximate conformal symmetry. 

The emission rate does not show the wiggly pattern (see Fig. 4) because the Planck factor 
usually suppresses the contribuation of higher partial wave except S-wave. The numerical 
calculation indicates that the emission rate generally decreases with increase of mass. 

Unruh has shown in his paper [16] that the S-wave cross section of Dirac fermion is exactly 
1/8 of that of scalar in the low-energy limit. It is interesting to apply our numerical method 
to check whether or not this ratio is maintained in the full range of energy. In Ref. [20] we 
examined the effect of mass in the absorption problem of the higher-dimensional object such 
as p-brane. Since, however, the low-energy of the scalar field is assumed from the beginning 
for the analytic computation, it was hard to understand the high-energy behavior. Recently, 
the absorption and emission problems of higher-dimensional black hole for massless particles 
with spin was numerically investigated in Ref. [34,35] in the full range of energy. In this 
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context it seems to be interesting to apply our method for the investigation of mass effect 
in the problem of higher- dimensional objects. 

Acknowledgement: This work was supported by the Korea Research Foundation Grant 
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Appendix 

In this Appendix we will show Eq.(3.5) explicitly. Since J~t(±\{ solutions, they 

should satisfy Eq.(3.1) separately. Multiplying Ft(±) to the equations of Tt^) and subtract- 
ing them subsequently, one can derive easily 

ax x[x — x s ) 

where W is the Wronskian defined in Eq.(3.5). Solving (A.l) easily gives 

W = C (A.2) 
x[x — x s ) 

where C is an integration constant. 

In order to compute C we use the explicit expressions of Tt{±) 

Fi i+) (x,x.) =i *-i e -*[*+¥H*-*.l] Q + (A.3) 

Then the explicit computation of W with these expressions yields the following leading term 

W ~ 4 ~ ? ^ \ (A.4) 
ar x(x — x s ) 

where the last approximation comes from the fact that J r e(±) is an asymptotic solutions. 
Thus we have C — 2i. 

Same procedure can be applied to Q^(x,x s ) and Ge(x,x s ), which yields C = —2i\gz\ 2 x 2 s /v. 
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FIGURES 

FIG. 1. Plot of the partial absorption cross sections for massless scalar case. The black lines 
represent the partial absorption cross sections computed by analytic solutions and their analytic 
continuations. The red lines are result of Ref.[16]. The fact ao = 47r at k = indicates the 
universality of the low-energy absorption cross section for S-wave. The peak points of a? are given 
in Table I. 

FIG. 2. The ^-dependence of the effective potential (2.4) for massless case. The potential 
makes a barrier, which separates the asymptotic and near-horizon regions. The barrier heights are 
given in Table I. 

FIG. 3. The energy-dependence of the total absorption cross sections whenm = 0, 0.15, 0.3, 
and 0.5. The wiggly behavior indicates that each partial absorption cross section has a peak at 
different energy scale. Regardless of mass the total absorption cross sections approach to 277r/4 at 
high-energy limit, which reflects the approximate conformal symmetry. 

FIG. 4. The energy-dependence of the total emission rate when m = 0, 0.15, 0.3, and 0.5. The 
disappearance of the wiggly pattern indicates that the rapidly decreasing Planck factor suppresses 
the contribution of the higher partial waves. The scalar mass generally reduces the emission rate. 

FIG. 5. The r* -dependence of the effective potential (3.2) for t = (Fig. 5a), t = 1 (Fig. 5b) 
and i = 2 (Fig. 5c) . The scalar mass generally lowers the barrier height and makes the near-horizon 
regime more stable than the asymptotic regime. When m is larger than the critical mass m c , the 
potential peak becomes negative and cannot play a role of barrier any more. 

FIG. 6. The energy-dependence of the imaginary part of the phase shift for I = (Fig. 6a), 
1 = 1 (Fig. 6b) and t = 2 (Fig. 6c). As in the case of massless scalar j3g is a monotonically 
increasing function. However, the scalar mass generally decreases the increasing rate of fy. 
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FIG. 7. The energy-dependence of the partial absorption cross section for i = (Fig. 7a), 
1 = 1 (Fig. 7b) and 1 = 2 (Fig. 7c). The scalar mass makes a local minimum in ao when m is 
less than the critical mass. If m > m c , the partial absorption cross sections become monotonically 
decreasing functions regardless of I. 

FIG. 8. The mass-dependence of the low-energy absorption cross section for S-wave. The 
increasing behavior of ao indicates that the universality is generally broken when the scalar field 
has a mass. The red line is a result of Ref.[16]. 
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Fig. 6b 
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Fig. 6c 
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